White noise flashing Brownian pump 
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A Brownian pump of particles powered by a stochastic flashing ratchet mechanism is studied. 
The pumping device is embedded in a finite region and bounded by particle reservoirs. In the 
steady state, we exactly calculate the spatial density profile, the concentration ratio between both 
reservoirs and the particle fiux. A simple numerical scheme is presented allowing for the consistent 
evaluation of all such observable quantities. 

PACS numbers: 05.40.-a, 05.70.Ln. 



I. INTRODUCTION 

Free diffusion by itself is not the appropriate phys- 
ical mechanism for selective transport processes. The 
appearance of net directed motion needs the breaking 
of detailed balance and spatial inversion symmetry. In- 
spired on Feynman's ratchet and pawl device such 
phenomenon has been named after the ratchet effect. 
Moreover, when the environment is thermally fluctuat- 
ing, Brownian transport is crucially affected and con- 
trolled by noise (for a broad review see and references 
therein). Its relevance to molecular engines, such as mo- 
tors, pumps and channels, has been proposed from the 
foundation idea of the so-called Brownian motors 3] . 

Brownian pumping is then an active nonequilibrium 
transport process in which fluctuations play a very im- 
portant role. Regarding experiments, Na,K-ATPase 
pumps have been perturbed by an oscillating electric field 
[3, [H| , driving ions whose net flux was measured as a func- 
tion of the amplitude and frequency of the field. The 
flux of particles created by pumping machines has been 
studied theoretically in @, 0, S @1- A typical ratchet 
mechanism is often assumed in the modeling and so nor- 
malization of the probability distribution and periodic 
boundary conditions are imposed. In this work, the scope 
is different. While the ratchet mechanism is still the re- 
sponsible for the transport of particles, we focus on the 
nonequilibrium concentration gradient created and main- 
tained by the pump at both sides of the membrane. The 
system is considered finite and not infinitely periodic. 
Then, our main observable is not only the flux J of par- 
ticles but the density profile p{x) in the membrane and 
the ratio of concentrations pi / po created at the particle 
reservoirs. See figure [1] 

The structure of paper is the following. First we in- 
troduce the model for the Brownian pump by means of 
a Langevin equation, which can be mapped into a See- 
beck ratchet |10l . [Til . [T3 |. In the steady state, a stan- 
dard theoretical analysis is carried out to obtain the non- 
equilibrium density of particles p{x) that is generated at 
zero and finite fiux. The flux J itself is also determined. 
Then we present a simple numerical scheme from which 
the densities and fluxes can be measured and satisfac- 
torily compared to the predictions. For a piece-linear 



Figure 1: Scheme of the pumping device: a spatially asym- 
metric and time dependent potential V{x,t), embedded in a 
finite region of length L (denoted as the membrane) flashes 
in time creating a density profile between the reservoirs of 
densities po and pi. 



saw-tooth potential we calculate explicitly the expres- 
sions derived in the analysis and explore them as a func- 
tion of the main parameters of the system. We end with 
some conclusions and comments for future work. 



II. GENERAL THEORETICAL ANALYSIS 

We consider an underdamped Brownian particle mov- 
ing under a time dependent potential V{x,t). The cor- 
responding equation of motion for its position is the 
Langevin equation 



= --fX-V'{x,t)+7j{t), 



(1) 



where m is the mass of the particle, 7 introduces the fric- 
tion and ri{t) is the thermal noise accounting for thermal 
fluctuations of the environment with the usual autocor- 
relation 



< f]{t)v{t') 2jkBTS{t - t'). 



(2) 
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The time dependent potential V{x, t) consists of a ratchet 
part V{x) (spatially asymmetric potential) time modu- 
lated by a stochastic process in the following form, 



V{x,t)^V{x)[l + x{t)], 



(3) 



where xit) is another white and Gaussian noise with zero 
mean, uncorrelated with ri{t), whose autocorrelation is 



<xit)xit') >^2Q5{t-t'). 



(4) 



In the regime in which friction dominates inertia, the 
Langevin equation reduces to its overdamped limit. 



x^-V'{x)^V'{x)x{t)+v{t). 



(5) 



Without loss of generality, the friction constant 7 has 
been absorbed in the time units. The last equation has 
a multiplicative noise and it should be treated with care. 
One must dilucidate first whether it has to be interpreted 
according to either Ito or Stratonovich rules. The order 
of the limiting procedures from which one arrives to an 
overdamped equation from an underdamped equation, 
having assumed first that the noise x(t) is obtained as 
the limit of an Orstein-Ulhcnbeck process, determines 
that th e appro priate stochastic interpretation is that of 
Ito [Tsl . Il4l . fisj . Then, equation (O can be rewritten as 



where 



x = -V'{x)+g{x)£,{t). 



l{x) = y/keT + QV'{x)\ 



(6) 



(7) 



and the new effective white noise ^ {t) has zero mean and 
correlation < £,{t)£,{t') 2S{t - t'). The continuity 
equation for the density of particles p{x,t) (the corre- 
sponding Fokker-Planck equation) is 



dtp{x,t) = -d^J{x,t), 



(8) 



which, by using Ito's prescription in ([S]), gives the follow- 
ing explicit expression for the flux J{x,t) 

- J(x, t) = V'{x)p{x, t) + [g{xfp{x, t)] . (9) 

In the steady state, the density is just a function of 
space and thus the flux becomes a constant, J. The 
density p{x) follows a first order non homogeneous linear 
differential equation, whose formal solution is 



p{x) = Z(x, Co, J)exp 
with 

Z{x,Co,J) = Co- J 



g'^{z) g{z) 



dz 



dz 



-2Sll4l)dz 

s(« ) / 



(10) 



(11) 



The unknown constant cq is found by imposing the left 
reservoir concentration po = p{0~), as a fixed bound- 
ary condition. Then cq = po- At each of the membrane 



boundaries with the reservoirs, this is at a; = and x = L, 
we will distinguish between approaching from left and 
right side because of the possible existence of discontinu- 
ities, which we will see in short. In what follows, we will 
study two different situations. 

First, we impose a zero total flux: J ~ 0. This cor- 
responds to the case in which the pump in maintaining 
the maximum concentration difference between the two 
reservoirs across the membrane with no net leaking of 
particles. This situation is analogous the stalling force in 
Brownian motors. From (jlOp . the density profile in the 
membrane is 



p{x) = p{xo) 



exp 



" V'jz) 
9Hz) 



dz 



(12) 



The exponent 2 in the prefactor before the exponential 
is a characteristic of the Ito's interpretation. It changes 
to 1 for Stratonovich's. One can find the expression for 
the ratio of concentrations at both sides of the pumping 
mechanism, being pi = p{L^). Assuming no system- 
atic drift in the system, this is V'{L~^) = V'{Q~), then 
g{L+) = g{0~) and therefore yields 



Pi 
Po 



— = exp 



1 



knT 



-V'iz) 



i + ^v'{zr 



-dz 



(13) 



We will study such expressions in more detail in subse- 
quent sections. 

The concentration po was fixed in the former case. If 
we now impose an arbitrary pi too, then the flux J cannot 
longer be zero. From (fTO|) and (fTTj) . one can get the 
explicit expression 



Po - Pie 



v'(^) I n a'(^) ' 



dz 



lo 



L + 



r (- 



(z')^ 



(14) 



dz 



Substituting J back in pT|) we can obtain, from ^TU\\ . the 
density p{x) in the steady profile for any choice of po and 
Pi- 

Finally, let us stress that although the resolution of 
the Fokker-Planck equation is standard, the conditions 
we have imposed are not the usual ones in ratchet models, 
leading to new solutions. The system is not infinitely pe- 
riodic and the probability distribution function is viewed 
as a density of particles and thus, it is not normalized. 
The present scheme is a realistic one for modeling pumps 
and channels operating between particle reservoirs. 

III. NUMERICAL SCHEME 

The time discretization and simulation of the particle 
dynamic evolution equations, ([5]) or ([21), are standard. 
The simplest method consists in using Euler's algorithm 
so that, for instance, from equation ([5]), the evolution of 
the position Xi{t) of the i~th particle turns out to be 

x,{t + 6t) = x,{t) - V'{x,{t))dt + g{x,{t))X,{t), (15) 
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Figure 2: Density profile p(x) in the steady state at J = 
obtained from numerical simulations. We use the potential 
introduced in Section IV. Reflecting boundary conditions are 
implemented at the end of both reservoirs. 

where 5t is the (small) integration time step and the 
stochastic term Xi (t) is constructed as, 

rt+St 

Xi{t)^J £,{t')dt' ^ V2Si a,, (16) 

in which tti are gaussian random numbers N{0,1). In 
fact, algorithm (jlSp corresponds to Ito's interpretation. 

The implementation of the boundary conditions needs 
a more careful study. The membrane is surrounded by 
the reservoirs, which we also include virtually in the sim- 
ulations, by means of an artificial length Iq. Then, the 
spatial domain of the system consists of three regions: 
the left reservoir x e [—Iq, 0), the membrane x e (0, L) and 
the right reservoir xe{L,L + Iq]. We consider a large 
number N of non- interacting particles, each one evolving 
according to (|15|) . When a particle is in a reservoir then 
V{x) = and it diffuses freely; only random kicks due to 
thermal noise drive the particle. Each particle is num- 
bered and its position, Xi{t), recorded. The numerical 
simulation always starts with particles distributed ho- 
mogeneously. The measures are taken waiting a certain 
time so that the system has reached the steady state. In 
general we observe a transient during which po(t) de- 
creases and pi(t) increases up to the steady state, in 
which pi > pq and it is kept constant. To obtain the 
density profiles p{x) an histogram of the distribution of 
particles is drawn. 

For the physical case of zero flux {J — 0), we impose 
reflecting boundary conditions at the points x — —Iq and 
a; = L + ^0 of our simulation system. For non zero flux 
(J 7^ 0), we take periodic boundary conditions at points 
— ^0 and L + Iq, ensuring the same flux along all the re- 
gions. These two simple rules lead to the physical sit- 
uations of interest, as we will show now. In the former 



Figure 3: Density profile p{x) in the steady state at J 7^ 
obtained from numerical simulations as in previous figure. 
Now periodic boundary conditions are applied at the ends of 
the reservoirs leading naturally to a constant J and fixed pi 
and po. 

case (reflecting boundary conditions), the steady state 
density of particles at both reservoirs is constant, as ex- 
pected for J = 0. See Fig. [21 One can plot the detailed 
shape of p{x) in the membrane and also measure the con- 
centration ratio pi/po that the pumping mechanism has 
created. We will analyze such observable quantities in 
detail in the following sections. The difference of con- 
centrations created and maintained implies that the flux 
was not zero during the transient and that the system is 
now in a nonequilibrium steady state. 

When imposing periodic boundary conditions in the 
simulations (not in the real physical Brownian pump 
since po ^ pi) we observe, in the steady state, a linear 
density profile in the reservoirs, signature of the finite flux 
J. Such flux should be constant and equal everywhere. 
In Fig. [3] a typical histogram is shown. This simple 
and intuitive method is able to generate data of J as a 
function of the boundary conditions; the concentrations 
Po and pi. More advanced and sophisticated studies of 
simulation of Langevin trajectories with specific bound- 
ary conditions can be found in Refs.[l3, H, EH. In our 
method, the pump is moving particles from the left to the 
right forcing a gradient of concentration in the reservoirs 
(which are connected only in the numerical scheme). Ac- 
cordingly, a net flux appears, which fulfills Fick's law so 
that 

J='^(.Pi-po)- (17) 

Therefore, independently of the initial conditions the 
system will evolve to the steady state in which the value 
of J in equations (fH)) and pT)) coincides. Then the flux 
can be measured by counting the net number of parti- 
cles crossing the point x = L + Iq (where the periodic 
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condition is implemented) as a function of time. A com- 
plementary way is to fit the histograms in the region 
(L, L + lo) U(^^o, 0) with the linearly decreasing behav- 
ior predicted by Pick's law. From the slope, J is found 
and also the concentrations pq and pi can be determined. 
The agreement between theory prediction for J{po,pi) 
and simulation results is good as we will show in the 
next section. 



IV. ANALYTICAL AND SIMULATION 
RESULTS 

In this section we complete the analytical results of 
Ref. [2Q| and compare them with numerical simulations. 
The explicit model we consider is a piece-linear saw-tooth 
potential depicted in Fig. [TJ It is defined in two regions 
as 



VAix) = Vo— xeiO,\L), 



'XL 



Vb{x) - Vo 



L 



(1 - \)L 



X e {\L, L). 



(18) 



(19) 



Vq is the height of the potential, A (which can only take 
values between and 1) controls the asymmetry and L is 
the total length where the pumping device is allocated. 

The forthcoming subsections are devoted to calculate 
and discuss the exact analytical expressions of the pro- 
file density of particles p{x) and the concentration ratio 
Pi/po at J = 0, as well as the normalized flux J/po as a 
function of pi/po- Since the potential is linear in pieces, 
the force is discontinuous, yielding to discontinuities in 
the density profiles, which we will carefully characterize 
without any further problems. 



A. Density profile p{x) at J = 

Let us recall expression for the spatial density pro- 
file p{x) given a baseline concentration po and zero flux 
conditions, J = 0. From the potential in (fT5|) and 
p{x) can be obtained exactly. Nevertheless as the poten- 
tial is linear in pieces, we expect discontinuities in p{x). 
Thus, we have to evaluate the profile in four different 
zones denoted by the subscripts 0, A, B and 1, whose 
meaning is clear from Figs. [1] and [H First we intro- 
duce the following dimensionless parameters to simplify 
calculations. 



Vo 



Vo 



QksT 
L2 ' 



(20) 



where vq is the relative energy barrier of the potential 
compared to thermal energy. The parameter a is measure 
of the strength of the flashing mechanism. 

The density of particles in region A, this is pA^x), in 
which X e {0~^,XL~), is obtained from (|12p by taking 




Figure 4: Density profile p{x) for no = 8, A = 1/3 and a = 
0.01 in logarithmic scale. The histogram from simulations 
(in gray) falls just on top of the theoretical predictions (solid 
lines) . 



xq — 0~ , so that p{xo) = po and gixo) = ^/ksT. Notic- 
ing that g{x) in region A (denoted by gA{x)) is a con- 
stant, one has 

gA = .g(0+) = g{XL-) = ^ ksT + Q {Vol Xhf (21) 
which finally yields to 

(22) 



Pa(,x) = — exp — 



where, for simplicity in the notation, the following new 
dimensionless quantities have been defined. 



LJl=l+a(y) , LJ2 = 1 



Vo 



(23) 



1- A 

Then, the jump of the concentration at a; = is simply 

1 



<5i=p(0-)-p(0+) = po 1- 



(24) 



The density of particles in region B is found similarly. 
We recall (IT^ and now choose xq = . From the 
above expressions and by noticing again that gsix) is a 
constant. 



gs = keT + Q {Vo/il ~ \)L f, 

the concentration profile Pb{x) is found to be 

Po f vo\ f Vo X- XL 
Pb[x) — exp exp — — - 

0J2 \ OJl J \ijJ2 L - XL 

The jump between zones A and B (at x — XL) is 
52 = piXL+) - p{XL-) =. pb{XL) - pa{XL) 



(25) 



(26) 











\W2 





(27) 
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Figure 5: Logarithm of the ratio of concentrations versus the 
asymmetry parameter A at J = 0. Note that it is antisym- 
metric under the transformation A ^ 1 — A. The parameters 
are vo = S and a — 0.01. Sohd line corresponds to theory, 
equation (|28p . and circles to simulation data. 



Figure 6: Ratio of concentrations pi/po as a function of the 
relative potential barrier «o at J = 0. The values of the 
parameters are A = 1/3 and a = 0.01. Line and circles as in 
previous figure. 



The constant density pi at the other side of the mem- 
brane is discussed in detail in the next section. The jump 
S3, which corresponds to x = L, is obtained likewise from 
the difference between pi and pb{L). In Fig. |4]we show, 
in logarithmic scale, the above analytical predictions for 
p{x) calculated in every region and the corresponding 
jumps. Note that every piece of prediction fits perfectly 
to the histogram (in gray) built from the data. 

B. Ratio pi/po at J = 0. 

We focus now on the value of the ratio of concentra- 
tions at both ends of the membrane that the pumping 
Brownian device is able to create and maintain at the 
stalling regime. For the linear saw-tooth potential, equa- 
tion (fT^ gives 




Let us explore this result on the parameters vq, a and 
A. In Fig. Owe check the symmetry properties of our pre- 
diction with respect the parameter A. This figure shows 
the right-left inversion symmetry of the problem when 
we change A for 1 — A. Note that at A = 0.5, although 
there is a time modulation of the potential, the device 
does not pump because the spatial inversion symmetry 
is not broken and so there is no preferred direction. The 
more asymmetric the potential is, the greater the pump- 
ing capacity is achieved. This does not mean higher effi- 
ciencies with respect to energy consumption. This issue 
is not studied in this work. 



In Fig. [6]the density ratio is studied versus the relative 
energetic barrier vq. There is an optimum value which 
gives the largest difference. If the barrier height is small 
with respect the thermal energy, diffusive loses through 
the pump are important and the ratio decreases. Note 
that for Vq there is no concentration difference; pi = 
Pq. On the other hand, for very large values of vq, few 
particles can cross the barrier and get to the other side 
despite the flashing of the potential. Then the pumping 
decreases again. 

In Fig. [7] we show the ratio of concentrations versus 
a (we have varied Q in the a exploration). For low val- 
ues of Q, the potential barely changes in time, diffusion 
dominates and thus particles can scarcely be pumped. 
For a very strong flashing the potential is so often dis- 
torted that particles do not have time to cross through 
the membrane and they are again poorly pumped. In 
between both regimes, there is a region that enhances 
transport. Such optimal regime indicates that the flash- 
ing intensity Q can be tuned to be optimal. This is a 
common feature of flashing ratchets The maximum 
appeared in Figs. [S] and [7] can be compared qualitatively 
with the experimental results of A, ,5] for the amplitude 
and frequency of the flashing perturbation. 

Let us analyze in more detail the case a = 0, which 
gives po = Pi- This limit is physically interesting because 
it corresponds to QkBT/L^ — > 0. When the intensity in 
the multiplicative noise vanishes {Q = 0), the breaking 
of detailed balance does not occur and, so, the ratchet 
effect cannot take place. Moreover, another way to make 
a vanish is setting T — 0. We have to be very careful 
because T also appears in vq = Vo/A;bT. In fact, in 
the absence of thermal fluctuations, the flashing ratchet 
mechanism still works because the multiplicative noise 
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Figure 7; Ratio of concentrations versus parameter q at J = 
for vo = 8 and A = 1/3. Line and circles as in previous figure. 
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Figure 8: Flux J/po as a function of pi/po- The system pa- 
rameters are A = 1/3, vq — 8, a = 0.02 and L = 1. Solid line 
corresponds to theory, equation (|30|l . and circles to simulation 
data 



does all the job (breaks detailed balance and supplies 
fluctuations). Therefore, in the a exploration we have 
kept T different from zero. 



C. Flux J/po versus ratio pi/po 

We can recall expression for the total flux of par- 
ticles and rewrite it as 



J 



1 - Pi/ Pi 



stall 



Po _J_ 



(29) 



where p^*"" is the concentration in the right reservoir 
when J — 0, which is calculated in For the linear 

saw-tooth potential, the above formal expression can be 
explicitly expressed as 



kBT\ poe-""/" 



Pie 



-Vo/uJ2 



Ml +N2 



where 



AAi^^^Afl-e-^^/^O 
Vo \ J 



n.^'^ii-x) (1 



(30) 



(31) 



(32) 



Expression (|30p fulflUs the symmetry J ^ — J when A 
1 — A and po <-> pi, which means that if we take the 
mirror image of the set up, we should see the same flux 
going to the opposite direction. Note that when J = 
only the ratio of densities is relevant, while for J 7^ 0, 
both values are needed separately. There is a decreasing 
linear behavior of J with respect to pi , as it is clear from 



simple inspection of ([29]) and If pi > pf the flux 
is reversed because entropic forces surmount the pumping 
driving. 

In Fig. [8] the normahzed flux J/po is plotted as a func- 
tion of the density ratio pi/ po- The simulation points are 
successfully obtained from the method explained in the 
numerical scheme section. Increasing the length of the 
left and right reservoirs the value of the flux J decreases. 
We can measure it from the slope of the linear proflle of 
the density in such regions and also extract po and pi. 
Despite some small deviations due to the errors in the 
process of data obtention because of flnite statistics and 
non zero (5i, we can say that theory and simulations are 
in agreement. This conflrms the validity of the numeri- 
cal scheme proposed, which allows to measure flirxes and 
concentrations from simple numerical simulations which 
have a clear physical interpretation. 



V. CONCLUSIONS 

Brownian pumps, unlike typical Brownian motors, do 
not aim to create a net flux of particles in periodic bound- 
ary conditions, but instead to achieve and actively keep a 
density gradient between two reservoirs. We have studied 
a simple model in which the ratchet effect together with 
appropriate boundary conditions leads to such mecha- 
nism. It has been thoroughly characterized from the 
theoretical point of view with analytical exact results. 
Moreover, a new and simple numerical scheme has been 
proposed to measure concentrations and fluxes, faithfully 
reproducing all the theoretical predictions. The efficiency 
of such devices is not easy to investigate since one should 
analyze how much energy the fluctuating potential is in- 
serting into de pump and what is the energetic profit 
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taken out from such input. This issue together with the 
use of physical parameters in the biological scale is un- 
der study. The present work is then a starting point for 
modeling nanometric bio- machines, such as channels and 
pumps, which control the flux of particles across the cell 
membrane. 
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